Born expansion of the Casimir-Polder interaction of a ground-state 
atom with dielectric bodies 

un : 

Stefan Yoshi Buhmann*, Dirk-Gunnar Welsch 

o 

' Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universitat Jena, Max-Wien-Platz 1, 07743 Jena, Germany, Phone: +49- 
\ 3641-947100, Fax: +49-3641-947102 

o ■ 

^3 ! Received: date / Revised version: date 



> 
00 

r- 

o 
in 
o 

^ : 
^— > ■ 

G ■ 

3 : 



Abstract Within leading-order perturbation theory, 
the Casimir-Polder potential of a ground-state atom pla- 
ced within an arbitrary arrangement of dispersing and 
absorbing linear bodies can be expressed in terms of the 
polarizability of the atom and the scattering Green ten- 
sor of the body-assisted electromagnetic field. Based on 
a Born series of the Green tensor, a systematic expansion 
of the Casimir-Polder potential in powers of the electric 
susceptibilities of the bodies is presented. The Born ex- 
pansion is used to show how and under which conditions 
the Casimir-Polder force can be related to microscopic 
many-atom van der Waals forces, for which general ex- 
pressions are presented. As an application, the Casimir- 
Polder potentials of an atom near a dielectric ring and 
an inhomogeneous dielectric half space are studied and 
explicit expressions are presented that are valid up to 
second order in the susceptibility. 

PACS 12.20.-m Quantum electrodynamics - 34.50.Dy 
Interactions of atoms and molecules with surfaces; pho- 
ton and electron emission; neutralization of ions - 34.20.- 
b Interatomic and intermolecular potentials and forces, 
potential energy surfaces for collisions - 42.50.Nn Quan- 
tum optical phenomena in absorbing, dispersive and con- 
ducting media 



1 Introduction 

The forces of electromagnetic origin that arise between 
electrically neutral, unpolarized but polarizable objects 
are commonly known as dispersion forces PI2E1E1E]- 
They were first addressed within the context of quan- 
tum electrodynamics (QED) by Casimir and Polder [SJ 
[7] , who showed that they may be attributed to the vac- 
uum fluctuations of the electromagnetic field. In accor- 
dance with the different nature of the interacting objects, 
one may distinguish between three types of dispersion 
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forces, namely the forces between atoms — in the follow- 
ing referred to as van der Waals (vdW) forces, the forces 
between atoms and macroscopic bodies — in the follow- 
ing referred to as Casimir-Polder (CP) forces, and the 
forces between macroscopic bodies — in the following re- 
ferred to as Casimir forces. 

Dispersion forces play a major role in the understand- 
ing of many phenomena, and they can be a useful or 
disturbing factor in modern applications. Apart from be- 
ing crucial for the understanding of many structures and 
processes in biochemistry [S], they are responsible for the 
remarkable climbing skills of some gecko and spider 
species |10| ; the construction of atomic- force microscopes 
is essentially based on dispersion forces while they 
are also responsible for the problem of sticking in nan- 
otechnology I n particular, CP forces between atoms 
and macroscopic bodies are needed for an understand- 
ing of the adsorption of atoms and molecules to surfaces 
|13| ; they can be used in atom optics to construct atomic 
mirrors |14j , while they have also been found to severely 
limit the lifetime of atoms stored on atom chips |15) . 

The study of CP forces which were first predicted 
for the idealized situation of a ground-state atom inter- 
acting with a perfectly conducting plate [Sj has since 
been greatly extended. Various planar geometries like 
the semi-infinite half space |16U17lll8lll9j , plates of finite 
thickness [201; two-layered plates [2J or planar cavities 
|20U22j have been treated, the most general planar ge- 
ometry being the planar multilayer system with an ar- 
bitrary number of layers |23U24j . Systems with spherical 
|25U26j or cylindrical symmetries |25U27| have also been 
considered. It should be mentioned that some theoretical 
approaches (in particular, those based on normal- mode 
quantization, e.g., Refs. |Btl8H9l2UI22l25p require a sep- 
arate treatment for each specific geometry, whereas oth- 
ers (in particular, the methods based on linear response 
theory, e.g., Refs. [TH1|17II21II26II27| ') lead to general ex- 
pressions that are geometry-independent. 

Recently, the problem has been studied within the 
frame of macroscopic QED in dispersing and absorbing 
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media and an exact derivation of a general expression 
for the CP force has been given 23 28 29 . Although the 
problem of calculating CP forces (or equivalently, the re- 
spective CP potentials) is thus formally solved, explicit 
evaluation requires knowledge of the (classical) Green 
tensor for the body-assisted electromagnetic field, which 
is (analytically) available only for a very limited class 
of geometries. In particular, inhomogeneous bodies or 
bodies of exotic shapes have not yet been treated. Nev- 
ertheless, it shall be demonstrated in this paper that the 
general solution — in combination with a Born expansion 
of the Green tensor — may serve as a starting point for 
the systematic study of a wide class of geometries. 

Furthermore, the Born expansion helps making gen- 
eral statements about two fundamental issues regarding 
CP forces. First, it may answer the question of whether 
and to what extent CP forces are additive. Second, it can 
be used to clarify the microscopic origin of CP forces. It 
is known that up to linear order in the electric suscep- 
tibility the force between an atom and a macroscopic 
body that may be regarded composed of atom-like con- 
stituents can be obtained by summation of two-atom 
(microscopic) vdW forces [5,30 ; an analogous relation 
between the Casimir force and CP forces can be estab- 
lished 30 31 32 33 34 35 36 . It is also known that pair- 
wise summation fails at higher order in the susceptibil- 
ity |37U38| . where many- atom interactions begin to play 
a role |39ll 40,41 ,42,43 ; in fact, it has been shown that 
an infinite series of many-atom interactions must be in- 
cluded in order to derive the CP force between an atom 
and a semi-infinite dielectric half space microscopically 

m 

The article is organized as follows. In Sec.Elthc Born 
expansion of the CP potential of an atom placed within 
an arbitrary arrangement of locally, linearly, and causal- 
ly responding isotropic dielectric bodies is given. The 
results are then used to elucidate the relation to micro- 
scopic descriptions (Sec. El, and to study some specific 
geometries (Sec.0. Finally, a summary is given in Sec. El 



2 Born expansion 

Consider a neutral, non-polar, ground-state atomic sys- 
tem A such as an atom or a molecule (briefly referred to 
as atom in the following) at position r A which is placed 
in a free-space region within an arbitrary arrangement 
of linear dielectric bodies. The system of bodies is char- 
acterized by the (relative) permittivity e(r,uj), which 
is a spatially varying, complex-valued function of fre- 
quency, with the Kramers-Kronig relations being satis- 
fied. Within leading-order perturbation theory, the CP 
force on the atom due to the presence of the bodies can 
be derived from the CP potential (see, e.g., Ref. |28|) 

U A (r A ) = ^ [°° duu 2 a A (iu)TrG^(r A ,r A ,iu) (1) 
2tt Jo 



according to 

F A (r A ) = -V A U A (r A ) 
(V A = V rA ). InEq. 0, 
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is the ground-state polarizability of the atom in low- 
est order of perturbation theory [uj^ q = (E^ ~ E^/h, 
(bare) atomic transition frequencies; = (O^d^ln^), 
clectric-dipole transition matrix elements of the atom], 
and G^ 1 ) (r, r', iu) is the scattering part of the classical 
Green tensor of the body-assisted electromagnetic field, 

G(r,r',w) = G (0) (r,r» + G (1) (r,r',w) (4) 

[G(°) (r, r', uj), vacuum part], which is the solution to the 
equation 



V x V x — - e(r, uj) 



G(r,r» = a(r-r')l (5) 



(I, unit tensor) together with the boundary condition 

G(r, r', ui) — > for |r — r'| — > oo. (6) 

Suppose now that 

e(r,cj) =£(r,a/) + x(r,w), (7) 

with the Green tensor G(r, r',w), which is the solution 
to 



V x V x 



■e(r,uj) 



G(r,r» = *(r-rOl, (8) 



being known. A (formal) solution to Eq. can then be 
given by the Born series 

G(r,r» = G(r,r',o,) 



fc=i 



2k 



II / d 3 *iX(sj,w) 



x G(r,si,w) • G(si,s 2) w) • • • G(s fc ,r',o;), (9) 
as can be easily verified using Eq. (JHJ together with 

k 



V x V x — - e(r, uj) 



E 

k=l 



21: 



n / ^i^j.w) 



x G(r,si,w) • G(si,s 2 ,w) • • • G(s fe ,r',w) 



(-) x(r,w)G(r,r',w). 



(10) 



Combining Eqs. 0,0, and 0, we find that the CP 
potential can be expanded as 



U A {r A ) = U A (r A ) + J2 AkU A {r A ), 



(11) 



fc=i 
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where 



U a (va) = 



2tt 



du u 2 a.A{iu) Tr (ta, i"A,iu) 



is the CP potential due to e(r, u>), and 
(-l) k hu 



(12) 



A k U A (r A ) 

c 00 
x / duu 
Jo 



2irc 2k 



II / d3s iX(s 3 , iu) 
.3=1 J 



x Tr[G(iv4, si, iu) ■ G(si,s 2 , iu)--- G(s k ,r A , iu)] (13) 

is the contribution to the potential that is of fcth order in 
x(r, oj). The Born expansion of the CP potential as given 
by Eqs. (fTT)) ~ l(l"3)l can be used to (systematically) calcu- 
late the potential in scenarios where a basic arrangement 
of bodies for which the Green tensor is known is (weakly) 
disturbed, e.g., by additional bodies or inhomogeneities 
such as surface roughness. 

Let us apply Eq. (|11|1 - (|13|1 to the case of arbitrar- 
ily shaped, weakly dielectric bodies, so that we may let 
e(r, u>) = l, and hence 



Ua(ta) =J2A k U A (r A ), 



(14) 



k=l 



where A k UA{^A) is given by Eq. JT3Jl, with 
G(r, r', iu) = Gv(r, r', iu) 



1 

-in 



(;)"*■ 



P 



= \(-) 2 6(p)\ + H v (r,r\iu) 

3 \u/ 

being the vacuum Green tensor (see, e.g., Ref. 
where 



H v (r,r',m) 



2 

C < 



Awu 2 p 3 



\—)p p 



a(x) = 1 + x + x 2 , 
b(x) = 3 + 3x + x 2 



(15) 
B31). 

(16) 

(17) 
(18) 



(p = r — r'; p=\p\; p = p/p). Combining Eqs. (|To^ and 
ljl"5|) [together with Eqs. l(TB|l - l(TS|) ]. one easily finds that 
to linear order in x the CP potential reads 



h f°° 

UA(rA) = A 1 U A (r A ) = - on , / dua A 



(iu) 



x ld* SX (sM MrA -t /C \ (19) 



where 



g 2 (x) = 2e~ 2x (3 + 6x + bx z + 2x A + a; 4 ). (20) 



2 , o™3 , „4\ 



In this approximation the CP force is simply a volume 
integral over attractive central forces, as is seen from 



52 (W c ) 



r 

7? 



4? r 



Y [6g 2 (ur/c) - (ur/c)g' 2 (ur/c)] 
2x (9 + 18a; + 16a: 2 + 8a; 3 

1 \ x—ur/c 



-3x 



(21) 



(r = |r|, f = r/r). 

In the retarded (long-distance) limit, i.e. 



(22) 



M 



where r_ = minjlr^ — s| : x( s ) 7^ 0} is the minimum dis- 
tance of the atom to any of the bodies, o/^ =min{cj^Q|n 
= 1, 2, . . .} is the lowest atomic transition frequency, and 
w^j is the lowest resonance frequency of the dielectric 
material, the exponential factor in g 2 (x) effectively lim- 
its the it-integral in Eq. (|19fl to a region where 

a A (iu) ~ a A (0), x(s, iu) ~ x(s, 0), (23) 

so Eq. (fCT?|) reduces to 

X(s,0) 



AMta) 



hca A {0) 1 d 3. 



327r 3 e 
23hca A {0) 
647r 3 e 



da; g 2 (a;) 



In the nonretarded (short-distance) limit, i.e., 



and/or r + <C 



(24) 



(25) 



M 



where r + = maxjlr^ — s| : x( s ) 7^ 0} is the maximum 
distance of the atom to any body part, uj\ = max{a^ |n 
= 1, 2, . . .} is the highest atomic transition frequency, and 
w^Jj is the highest resonance frequency of the dielectric 
material, the factors otA(iu) and x( s > iu) effectively limit 
the it-integral in Eq. (|19|l to a region where x — u\ya— s|/c 
<C 1 , so we may set 



32(2:) ~ g 2 (0) = 6, 



(26) 



resulting in 



3h 



167T 3 £q 



du a A (iu) 



d 3 s x(^4 



(27) 

The second-order contribution A 2 U~a(ya) can be sep- 
arated into a single-point term and a two-point correla- 
tion term, 

A 2 U A (v A ) = A\U a {ta) + A 2 2 U A (r A ), (28) 

as can be seen from Eqs. (11311 and 115fl for k = 2. The 
single-point term 

h 



A\U A (r A ) = 



967r 3 e 



du aA(iu) 



d*s X 2 (sM 92HrA ~:l /C \ (29) 
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which arises from the ^-function in Eq. (|15|l . differs from 
the first-order contribution AiU A {v A ) according to 



X(r,w) 



(30) 



hence its asymptotic retarded and nonretarded forms 
can be obtained by applying the replacement 1|3U|I to 
Eqs. J2U) and J2ZJ, respectively. 

The two-point correlation term is derived to be 



A 2 2 U A (r A ) 



1287r 4 e 



dua A (iu) / d 3 Six( s ii*u) 



"(tMtMt 



,-ii(a+/3+ 7 )/c 



'U(3- 

c 

fua\ /u(3\ /M7\ fua\fuj3\ /try 

\-) a y-n-)- a \-) h \-n- 



where 

g 3 (u, a,/3,7) = < 

-6 

-(t)"(t)»(?)+'(t)'(¥)"(?)(*A 
+o( ^ >( ^ )(M) > 

+0(t>(t)<^> 2 

-(.(=)i.(f)i(^)(a.3)0-7)(T«) : 

with the abbreviations 

- - i i - - a 

a = r ,\ — Si, a = at , ot = — 

a 
/3 



^ = 81-82, 0=101, /3: 
7 = s 2 -r A , 7=|7|, 7 : 



7 

7 



(33) 
(34) 
(35) 



having been introduced [recall Eqs. I|17|) and i|18[l ]. Note 
that the two-point contribution to the CP force, Eq. i|31|) . 
is a double spatial integral, the integrand of which can 
be attractive or repulsive, depending on the angles in the 
triangle formed by the vectors a, (3, and 7. 

In the retarded limit, where the inequalities (|22[1 
hold, the u-integral is again effectively limited to a re- 
gion where the approximations Ij23(l are valid, so Eq. (|3I|I 
reduces to 



(36) 



a 3 /3 3 7 3 
We introduce the notation 

a, = a l + f3 l + Y, i=l,2,3 (37) 
and perform the u-integral with the aid of the relation 



(38) 



Exploiting the triangle formula 
T = l-(a./3) 2 -(/3- 7 ) 2 -(7- ( i) 2 + 2(a^)(/3-7)(7.a) 







(39) 



[which is a trivial consequence of Eqs. H33fl - H35[) ] by 
adding the expression 



6T +—{[a 5 (P + 1 ) + /3 5 ( 7 + a) + 7 5 (a + 0)] 

+ 7[a 4 ((3 2 + 7 2 ) + /3 4 ( 7 2 + a 2 ) + 7 4 (a 2 + /3 2 )] 
+ 12(a 3 ^ 3 + /3 3 7 3 + 7 3 a 3 ) + 12a/3 7 (a 3 + /3 3 + 7 3 ) 
+ 52af3j[a(3{a + (3) + f3j{(3 + 7) + 7"(7 + a)] 
+ 138a 2 /5 2 7 2 } (40) 

to Eq. (|36fl . the result may be written in the form 
A 2 2 U A {r A ) = J d 3 six(si,0) J d 3 s 2 



X(s 2 ,0) 



(41) 



a 3 /3 3 7 3 (a+/3+7) 

+ / 2 (a,/3,7)(/3.7) 2 + / 2 (/3,7,a)(7-a) 2 
+ / 3 (a,/3,7)(a./3)(/3-7)(7.a)], 



where 



( 32 ) f x (a, /?, 7) 



2 2 

9 - 39^| + 22^| + 54^| - 65^2- -I- 20^| 



07 



/ 2 (a,/3,7) = 3 



CTi 

r«,2 



01 



0"l 



a 2 3a 2 (/3 + 7 ) 



+ 



4/3 7 (3a 2 - /3 7 ) 20a/3 2 7 



07 



(42) 



(43) 



/3(a,/3,7) = 

2 2 

- f - 39^| + 17^| + 72^| - 75^^ + 20^| (44) 
af a\ a\ o\ a\ 

[recall Eqs. (EHJhGS) as well as Eq. $7\ ]. 

In the nonretarded limit, where the inequalities (|25|l 
hold, the u-integral in Eq. (|31|l is effectively limited to a 
region where 

g 3 (u,OL,/3,j) ~ 53(0, a, 0, 7) 

= 3[l-3(d-/3)(/3- 7 )(7-a)] (45) 

[recall Eq. (|3*^|l : note that max{|si — S2 1 : x(si) 7^ 0, 
X(sa) 7^ 0} < 2max{|r A - s| : x(s) # 0}], so Eq. (JS) 
reduces to 



3ft 



dua^ra) / d six(si,zm) 



128^% Jo 
:i , . , l-3(a-/3)(/3- 7 )(7-a) 



x / d 3 s 2 x(s2,iu) 
[recall Eq. 



a 3 /3 3 7 3 



(46) 
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3 Relation to microscopic many-atom van der 
Waals forces 

In order to gain insight into the microscopic origin of 
the CP potential as given by Eq. (QJ, let us suppose 
that the susceptibility x( r > w ) is due to a collection of 
atoms of polarizability afi(w) and apply the well-known 
Clausius-Mosotti formula (see, e.g., Ref. 0H]) 



£ 1 n(r)a B {u>) 
1 - 1 n(r)a B (uj) 



(47) 



where n(r) is the number density of the (medium) atoms 
[n(r) = for r = r A ]. Since x(r,w) is the Fourier trans- 
form of a (linear) response function, it must satisfy the 
condition 



%(r, 0) > x(r, iu) > for u > 0, 



(48) 



which, with respect to Eq. (|47|l . implies that the inequal- 
ity 

±e' 1 n(r)a B (iu) < 1 (49) 

must hold. Substituting the susceptibility from Eq. (|47|) 
into the Born series of the CP potential U A (r A ) as given 
by Eqs. (|1 l|l — . taking into account that the Green 
tensor G(r, r',iu) can be decomposed as 



for'.tu) = l(-) 2 6{p)\ + H(r,r',iu) 
3 \u/ 



(50) 



[recall Eqs. © and l|15fl ]. and recalling the inequality 
<|49|) . it can be shown after some lengthy calculation 
that the Born series can be rewritten as an expansion 
of U A (i"a) in terms of many-atom interaction potentials 
U A B...B{rA,s 1 , . . . ,si) (see App.|I}, 



U A (r A ) = U A (r A ) 

1 ^ f 

+ Eyi II / dV*(^) 



1=1 



U A B...B(r A , si, . . . ,si) 



(51) 



where 

Uab...b(*i, ■ ■ - ,r/+i) 
(-l)'Vo +1 



{1 + Su)tt Jo 
x 5Tr[H(ri, r 2 , iu 



du u 2l+2 a A (iu)a l B (iu) 

H(ri +1 ,n,iu)]. (52) 



Here the symbol S introduces symmetrization with re- 
spect to ri, . . . , r; + i according to the rule 

5 1 &[R(ri,r 2> w).-.H(r i ,ri,w)] 

= Tr [ H ( r 7r (i)> r 7r(2),w) • • • H(r ff(j) ,r w(1) ,w)]. (53) 

The sum in Eq. (|53|l runs over the maximal number of 
j!/ [(2-S 2j )j] permutations it € P(j) £ P(J) [P(j) be- 
ing the permutation group of the numbers 1, . . . , j] that 



cannot be obtained from one another via (a) a cyclic per- 
mutation or (b) the reverse of a cyclic permutation (cf. 
App. EJ- The potential U A b...b(?i, ■ ■ ■ , >*i+i) is nothing 
but the (microscopic) vdW potential describing the mu- 
tual interaction of a (test) atom A at position ri and I 
(medium) atoms at different positions r 2 , . . . r; + i. 

The very general equation |JT]|. which follows from 
QED in causal media, gives the CP potential of an atom 
in the presence of macroscopic dielectric bodies in terms 
of the atomic polarizability and the (scattering) Green 
tensor of the body-assisted Maxwell field, with the bod- 
ies being characterized by a spatially varying dielectric 
susceptibility that is a complex function of frequency. 
Equation l|51|) clearly shows that when the susceptibil- 
ity is of Clausius-Mosotti type, i.e., Eq. I|47|l [together 
with the inequality l|49l) ] applies, then the CP potential 
is in fact the result of a superposition of all possible mi- 
croscopic many-atom vdW potentials between the atom 
under consideration and the atoms forming the bodies. 
Note that for the vacuum case, e(r,u>) = 1, a similar 
conclusion has been drawn by combining normal mode 
quantization with the Ewald-Oseen extinction theorem 
|38| . Moreover, from the derivation given in App. E| it 
can be seen that when Eqs. (|51[) and (|52(l [together with 
the inequality i|49|)] hold, then the susceptibility must 
necessarily have the form of Eq. (|47|l . 

From the above it is clear that in order to estab- 
lish the identity (|51|l to all orders in as, it is crucial to 
employ the exact relation l|47|l between macroscopic sus- 
ceptibility and microscopic atomic polarizability rather 
than its linearized version %(r, uj)—eQ 1 n(r)aB(^), which 
is is known to be sufficient for finding a correspondence 
between macroscopic and microscopic potentials to lin- 
ear order in \ (or t\ B . respectively) |51l3(JU31ll32ll33| . It 
should be pointed out that in more general cases where 
the susceptibility is not of the form l|47l) the result of ap- 
plying the Born expansion cannot be disentangled into 
spatial integrals over microscopic vdW potentials in the 
way given by Eq. (|51|l together with Eq. I|52|) . Obviously, 
the basic constituents of the bodies can no longer be ap- 
proximated by well localized atoms. 

According to Eq. Q , the expansion in Eq. I|51|) does 
not necessarily refer to all bodies. Hence from Eq. (|52() 
it follows that the many-atom vdW potential on an ar- 
bitrary dielectric background described by e(r, uj) reads 

U Al ... Aj (ri, . . .,Tj) 

(-ly^Vo r , 2? , ,. , 

= ———— — / duu J a Al (iu) ■ ■ -a A Aiu) 
{l + 5 2j )n Jo 

x 5Tr[H(r 1 ,r 2 ,m)---H(r j) r 1 ,m)], (54) 

the derivation being unique when requiring the vdW 
potentials to be fully symmetrized. In particular, for 
j = 2 [where j!/[(2— &2j)j] = 1, so the sum in the r.h.s 
of Eq. (|53|l contains only one term and symmetrization 
is not necessary], Eq. H54|) agrees with the result that 
can be found by calculating the change in the zero-point 
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energy of the system in leading-order perturbation the- 
ory |46| . In the simplest case of vacuum background, i.e., 
e(r,w) = 1, H in Eq. (O becomes H v [recall Eqs. l[TS|l- 
l|18|) ]. leading to agreement with earlier results |41U42| . 

Bearing in mind the general relation (|51() between 
the CP potential and many-atom vdW potentials, ex- 
plicit expressions for the two- and three-atom vdW po- 
tentials on vacuum background can easily be obtained 
from formulas given in Sec. H together with Eq. l|4"7|) . 
From Eqs. (jT§jl . J53J, and J57J one can infer the well- 
known result |H] 



E/As(ri,r 2 ) 



du g2(ur / c)otA(iu)aB{iu) (55) 



r = |ri — r2 1 ) , recall Eq. (|20|) . which reduces to 

23Hca A (0)a B {0) 



U A B{ri,r 2 ) 
in the retarded limit, 



647r 3 £gr 7 



0J- 



(56) 



(57) 



with r_ = r and cj_ = min{o; A , cj b }, and to 
3ft 



U AB (r 1 ,r 2 ) = - 167r3£ 2 r a 7o 
in the nonretarded limit, 

c 

r+ < , 

with r + = r and w + = max{w|,oig}. 
Similarly, Eq. (|31[1 implies that 

H 



duaA(iu)aB{iu) (58) 



(59) 



fZABc(ri,r2,r 3 ) 



dw 



64 7 r% 3 r 3 2 r 3 3 r 3 1 y Q 
x aA{iu)a B (iu)a c (iu)g 3 (u, r 12 ,r 2 3,r 3 i) (60) 

(r ii3 ■ =rj - Tj, r tj = \rij\ for i,j = 1,2,3), recall Eq. J32), 
in agreement with the result found in Refs. |40II41II42| . 
Equations lj4*Tl) and l|4*rJ)l show that Eq. (|6T7|l simplifies to 

TT hca A (0)a B {0)a c (0) 
U A Bc(ri,r 2 ,r 3 ) - 3 3 3 3 — — — — — ^ 

107T £ r 12 r 23 r 3ll r 12+ r 23 + ^3lJ 

x [/i (^12 ,^"23 ,^31 ) + /2 (^31 1 n2 ,r 23 )(r 12 -r 23 ) 2 

+ h (l2 , ^23, r 3i ) (f 23 • f 31 ) 2 + /2 (^23 , »*31 , r i2 ) (? 31 • f 12 ) 2 
+ /3 (»"12 , f 23 , r 3 l ) (f 12 ' f 23) (f 23 ' ? 31 ) (? 31 • f 12 )] (61) 

[fjj = rij/rij, recall Eq. (J^ and Eqs. l(4^ )l -l(il )l ] in the 
retarded limit [Eq. (|57|) with r_ = min{ri2, r 2 3, r 3 i} 
and u>- = min{w^ , , loq }] , and reduces to the Axilrod- 
Teller potential 

J7ABc(ri,r 2 ,r 3 ) 

_ 3ft[l - 3(fi2-f 23 )(f2 3 -f 3 i)(f 3 i-fi2)] 



647r%3r 3 2 r3 3 r| 1 

/>oo 

x / &ua A {iu)aB{iu)ac{iu) 
Jo 



(62) 




Fig. 1 An atom near a dielectric ring (schematic picture) 



in the nonretarded limit [Eq. I|59|) with r + = max{ri2, 
^2 3 ,^3i} and w + Emax{w^,4,u^}]. 



4 Application to specific geometries 

4-1 Dielectric ring 

Let us use the Born series given in Sec. |2] to calculate 
the (leading contributions to the) CP potential for some 
specific geometries and begin with a ground-state atom 
placed on the symmetry axis of a homogeneous dielec- 
tric ring of susceptibility x( w )> having radius Tq, (circu- 
lar) cross section ira 2 (where a<C r$), and volume V = 
2ir 2 roa 2 , the atom being separated from the center of the 
ring by a distance z A (Fig.^l. From Fig. we see that 
\r A — s| ~ \J z 2 ^ + Tq = p A for a <C tq, so an evaluation 
of the (trivial) volume integral in Eq. I|19fl results in the 
first-order CP potential 



A!U a {pa) 



HV 



32w 3 e p A J 



dua A (iu)x(iu)g 2 {up A /c), (63) 



which is attractive, as expected. In the retarded limit 
[Eq. J22J with r- = p A \ Eq. (jHSJ reduces to 

AiUa{pa) = 6^eop A (64) 

[cf. Eq. (J23J] , while in the nonretarded limit [Eq. (|2"5T 
with r+ = pa] one easily finds 



AiUa(pa) = - 



3hV 



l&ir 3 e p A J 



du a A (iu)x(iu) (65) 



[cf. Eq. H27fl ]. In both limiting cases the CP potential 
thus reduces to simple asymptotic power laws in p A , 
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where as usual the leading (inverse) power is increased 
by one when going from the nonretarded to the retarded 
limit. 

The first (single point) second-order correction term 
A\Ua{pa) can simply be obtained from Eqs. 1)63(1 ((65(1 
by means of the replacement 1(30(1 . while the calculation 
of the second (two-point) term A\Ua{pa) is a lot more 
difficult due to the factor |si — Sa|. We find (see Add. IB)) 

aIUa(pa) = (^^■°W0)x'(Q) (66) 



TT 3 e p 7 A 



in the retarded limit and 



A\U a {pj 



(0.08±0.03)fty 



n 3 e p A 



du a A(iu)x 2 \iu) 



(67) 

in the nonretarded limit. Recalling Eqs. ((28(1 and 1(30(1 . 
Eqs. ((64(1 and 1(66(1 imply that up to quadratic order in 
X we have 



Ua(pa) = - 



23hcVa A {0)x(0) 
64tt 3 £ p 7 a 



1 - (0.47±0.05)x(0) 



(68) 

in the retarded limit, while Eqs. I(65|) and 1(6 7|) show that 
in the nonretarded limit 



Ua(pa) = - 



3hV 



16tt 3 s p a 

x / du a A (iu)x(iu) 1- (0.77 ± 0.17)x(m) . (69) 
Jo L - 1 

The uncertainty in the magnitude of the contribution 
quadratic in x is due to the approximations made when 
calculating A\U(pa) (cf. App.[5J|. However, irrespective 
of these approximations, Eqs. ((68|) and ((69(1 show that 
the leading non-additive correction A\X) a(pa) to the lin- 
ear result AiUa(pa) does not change the powers in the 
asymptotic retarded and nonretarded distance laws (p A r 
and p~a i respectively), but merely modifies the constants 
of proportionality. A similar result has been found when 
studying a dielectric half space 



4-2 Many-body decomposition 

The explicit evaluation of multiple spatial integrals [be- 
ing the main difficulty when evaluating /^[/(r^)] can in 
fact be avoided in many cases by an appropriate decom- 
position of the body of interest, as shall be demonstrated 
in the following. To that end, let us decompose the body 
described by x(r,w) [recall Eq. JJJ] into smaller bodies 
numbered by n, so that 



where 



lv„(r) 



1 for r e V n , 
for r f V n . 



(70) 



(71) 



Substituting Eqs. lO and JHJ into Eq. (JT3J), and slight- 
ly rearranging the terms, we obtain 



A k U A (r A ) = J2 A kUA(rA), 



(72) 



where 



with 



A{U A {v A )= J2 AuT- ni ^A) (73) 



ll<....^rt; 



A l k U%- n '(r A ) 



(mi,...,rn fc )eX^ 1 



w ™ t ... mk 



(r A ) (74) 



is the sum of all /-body contributions of order k in x- I n 
Eq. O, 



(-i) k hpo 

2irc 2k 



< I dim + qa(w) 
'o 



II / d 3 SjX mj {sj,iu) 
j=i v mj 



x Tr[G(r A ,si, iu) ■ G(si,s 2 , iu) ■ ■ ■ G(s k ,r A , iu)], ( 75 ) 
and the notation 



-"ni...ni 



| (mi, . . . ,ra fc ) £ {m, ■ ■ ■ ,ni} k \ Vi 3j : mj = n; j 

(76) 



is used. 

In particular, to linear order in x we have 
A l U A {v A ) =A\U a {v a ) = Y,A\Ua 1 {*a) 

n 

= J2 W a(^a), (77) 

n 

so the CP potential is additive in this order. For the term 
quadratic in x, Eqs. l(72 ^ -l|71 )l together with Eq. (O 
reduce to 

A 2 U A (r A ) = AIU a {t a ) + A 2 2 U A (r A ) 

= j2^(T A ) + j2^ur(r A ) 

n m < n 

= J2 wt{va) + J2 [ w A n (*A) + W2 m (r A )] ■ 

n m<n 

(78) 

Obviously, the second term on the r.h.s. of Eq. 1(781) is 
the (overall) two-body contribution to the CP potential 
up to quadratic order in x- It can be regarded as being 
the leading correction to the additivity of the potential. 
Clearly, the expansion can in principle be extended to 
arbitrarily high orders in x, whereby fc-body interactions 
first appear at fcth order in x- I n particular, Eqs. ((77|) 
and (|75|) generalize the result that up to linear order 
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in x the CP potential of an atom near a homogeneous 
semi-infinite dielectric half space can be written as an 
(infinite) sum of thin-layer potentials |23| . whereas the 
contribution quadratic in x also contains two-layer terms 
leading to a breakdown of additivity |24|. 

To illustrate the application of Eqs. fT^ - lfTH)) to the 
calculation of the CP potential of complex bodies via 
decomposition into simpler bodies, let us consider an 
atom at position — za (za > 0) near a semi-infinite half 
space filled with a stratified dielectric medium, i.e. 



(79) 



where p(z) is some profile function [p(z) > for z > 0, 
p{z) = for z < 0], which may be normalized such that 
maxp(z) = 1. We decompose the half space into a set of 
thin slices of equal thickness d such that 



i{p'(z)\z > 0} < 1. 



(80) 



From Eq. l[77jl it follows that A\Ua(za) is given by the 
sum over the slices, which for d^C za contribute 



A\Ua(za) = ~~^~T I duu 2 a A (iu) J dqq 



x e 



-2b(zA_-\-na 



1 



If U 

2\Vc 



where 



X(iu)p(nd), (81) 



(82) 



(cf. Eq. (71) in Ref. PU), so after turning the sum into 
an integral one obtains 



A\U a {za 



hfio 



where 



47T 2 Jo 

dqqe- 2bZA P(2b) 



P(x) 



duu aA{iu)x{iu) 

2 



bcX 1 / u 
~u) ~ + 2 \6c 



2l 



dze xz p{z) 



, (83) 
(84) 



is the Laplace transform of the profile function p(z). In 
a similar way, the results 



AlUUzA) 



h/igd 
4tt 2 



duu 2 aA(iu) / dqq 



x e 



-2b(z A +nd) 





'bcX 


3 


4( 


< 'u 




V U J 


-4 H 




\ be 



2l 



X 2 (iu)p(nd) 

(85) 



(cf. Eq. (74) in Ref. H3j) and 
hu^d 2 



A 2 U% n 



(za) 



2tt 2 



duu 2 aA(iu) / dqqb 



x e 



-2b(z A +nd) 



"1 










2 











x x{i u )p{ m d)x{i u )p( n d) 



(86) 



(cf. Eq. (75) in Ref. [23]) can be derived, leading to 

hfj, 



A\U A {z A ) 



it dqqe 
and 



(i 



4tt 2 
- 2bZA P(2b) 



duu aA(iu)x (iu) 







3 








v U J 


-4 H 


4( 


6c 



A%U a {za 



hfio 

2tt 2 
dq q6e 

dzp(z) 



du u 2 aA{iu)x 2 {iu) 



(87) 



26z A 


'1 












2 








,6c J J 



dz'e~ 2f,z p(z') 



respectively. Hence the CP potential of the inhomoge- 
neous half space up to quadratic order of x has been 
calculated from the known CP potentials of one and two 
thin plates of constant permittivities. Needless to say 
that the method can be carried out to higher orders of 
X and can also be applied to other than planar systems. 

To give an example, we consider a dielectric medium 
whose permittivity oscillates in the z direction, 



e(r, w) = 1 + x( w ) cos 2 (fc z 2:)6'(z) 



(89) 



[&(z), unit step function]. Using Eqs. ffity. l|S7l) . and 
Q88[l. we find that up to quadratic order in x the CP 
potential takes the asymptotic form (see App. ICl 



Ua(za) 



~\A 1 C 4 F 3 [k z ZA) 



A2C4 



126 
169 



F 3 (k z z A ) + l %H 3 (k z ZA)]} (90) 



in the retarded limit [Eq. (1221 with r_ = za] and 
(A 1 C 3 + A 2 C 3 )F 2 {k z z A ) 



Ua(za) = — 

in the nonretarded limit [Eq. (|25|1 with 
23Hca A (0)x(0) 



(91) 



AtC A 
A 2 C 4 



AxC 3 



32tt 2 £ Jo 



6407T 2 £ 

169ftca A (0)x 2 (0) 
89607r 2 e 



dw ceA(iu)x(iu) 



za}- Here, 
(92) 
(93) 

(94) 



A 2 C 3 = - 



64tt 2 £ J 



duaA(iu)x 2 (iu) (95) 



are just the linear and and quadratic expansions of the 
well-known coefficients for the homogeneous half space 
in the retarded and nonretarded limits [where Ua(za) 
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-, 
3 




2.25 



Fig. 2 The CP potential of a ground-state atom in front of a 
half space with oscillating susceptibility in the retarded limit 
is displayed as a function of the distance za, where k z c/uj^ 
— > oo (upper solid line), k z c/u~^ —4 (dashed line), k z c/uo^ = 
2 (dash-dotted line), k z c/u^ = 1 (dotted line), k z c/io^ — > 
(lower solid line), x(0) = |. 



= — C±/z\, and = — C$/ z\, respectively], and the 

structure functions 



FM = ^ / 
]■ Jo 

H 3 (x) = - / 
J ■ Jo 



dtth 



dtVe 

are normalized such that 
H 3 {x) 



21 



2t 2 + x 2 
t 2 + x 2 ' 

2t 6 + 8x 2 t 4 + 5xH 2 + 2x 6 
(t 2 +x 2 ) 2 (t 2 + 4x 2 ) 



1 


for x — 




1 

2 


for x - 


-> oo, 


1 


for x - 




1 

4 


for x - 


-> oo 



(96) 
(97) 

(98) 
(99) 



Equations (|90|l and (|91|l . respectively, are illustrated 
in Figs. |21 and It is seen that the potential curves for 
different values of k z lie between the two solid curves that 
correspond to the limiting cases k z — > oo (upper curves) 
and k z — > (lower curves, which represent the potential 
observed in the case of the respective homogeneous half 
space). From Eqs. (|9"Uj) . lj$T)) . and JjJSJ it follows that the 
upper-curve potential values obtained in linear order of 
X are 1/2 times the lower-curve ones, which reflects the 
fact that for k z — > oo the potential in linear order of \ 
is simply determined by the average permittivity e(r, iS) 
~ 1 + \x{lo)0{z), cf. Eq. JBSJ- The factor found for the 
quadratic-order term is equal to 1/2 in the non-retarded 
limit [cf. Eqs. QTJ and but equal to 295/676 in the 

retarded limit [cf. Eqs. (O, lEHJ, an d lEU], owing to the 
influence of the two-plate term A^Ua^a)- Note that the 
curves for the intermediate values of k z approach the up- 
per limiting curve for large values of za and the lower 
limiting curve for small values of za, the potentials thus 
being (near za — k^ 1 ) somewhat steeper than z^ 4 and 



3 



< 

+ 




-100 



zaujX/c 



Fig. 3 The CP potential of a ground-state atom in front 
of a half space with oscillating susceptibility in the non- 
retarded limit is displayed as a function of the distance za, 
where k z cjuj\ — » oo (upper solid line), k z c/u)\ =20 (dashed 
line), k z c/u>\=& (dash-dotted line), k z c/ui\ = 2 (dotted line), 
k z c/u>2 ~*0 (lower solid line). 



z A — the power laws observed in the case of a homoge- 



neous half space. By controlling k z 
control the shape of the potentials. 



one can therefore 



5 Summary 

Within leading-order perturbation theory, the CP po- 
tential of a ground-state atom near dielectric bodies can 
be expressed in terms of the atomic polarizability and 
the scattering Green tensor of the body-assisted electro- 
magnetic field, where the bodies are characterized by a 
spatially varying dielectric susceptibility that is a com- 
plex function of frequency. Starting from this very gen- 
eral formula, we have performed a Born expansion of the 
Green tensor to obtain an expansion of the CP potential 
in powers of the electric susceptibility. The expansion 
shows that only in linear order the CP force is a sum 
of attractive central forces, while higher-order terms are 
unavoidably connected to multiple-point correlations in 
the dielectric matter, leading to a breakdown of additiv- 
ity 

Using the Born series, we have shown that when the 
dielectric bodies can be described by a susceptibility of 
Clausius-Mosotti type, i.e., when the basic constituents 
can be regarded as atom-like, then the CP potential is 
the (infinite) sum of all microscopic many-atom vdW 
potentials between the atom under consideration and the 
atoms forming the bodies. As a by-product, a general 
formula for the many-atom vdW potential of arbitrary 
order and on an arbitrary background of dielectric bodies 
has been found, which generalizes previous results found 
for atoms in vacuum. 
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Apart from being useful for making contact with mi- 
croscopic descriptions of the CP force, the Born series 
can also be used for practical calculations, particularly 
when the Green tensor is not available in closed form. We 
have employed two strategies, (i) By direct evaluation of 
multiple spatial integrals, we have determined the at- 
tractive CP potential of a weakly dielectric ring, finding 
asymptotic l/p A and l/p A power laws in the retarded 
and the nonretarded limits, respectively, (ii) By reduc- 
tion to simpler bodies with known CP potentials, we 
have derived expressions for the CP potential of an atom 
placed in front of an inhomogeneous stratified half space, 
with special emphasis on an oscillating susceptibility. In 
this case the potential exhibits — for distances compara- 
ble to the oscillation period — a somewhat stronger power 
law than in the case of a homogeneous half space. 
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A Derivation of the expansion in terms of 
many-atom vdW potentials [Eq. (15 111 ] 

As a preparation, we derive the symmetrization 1531) . 
The completely symmetrized form of a many-atom po- 
tential is given by 

S/(n, ....I-,-) = 1 /(Mi)>-- •>**(,•}), 

V 23 U TreF(j) 

(100) 

where P(j) denotes the permutation group of the num- 
bers 1, . . . , j and 1/(2 — 52j)j is a normalization factor. 
As a trivial consequence of the cyclic property of the 
trace as well as the symmetry property of the Green 
tensor 03] 

G(r,r» = G T (r',r,c) (101) 
together with Eq. 1)50(1 one easily finds that 

Tt[H(ri ) r 2 ,w)---H(r J ,n,w)] 

= Tr[R(r 7r ( 1 ),r 7r( 2),a;) ■ • • H(r„^,r„^,u)] (102) 

if 7r is either a cyclic permutation [e.g., 7r(l) = 2, 7r(2) = 3, 
. . . , 7r(j) = 1] or the reverse of a cyclic permutation [e.g., 
7T(1) = j, 7T(2) =j — 1, ... , 7T(j) = 1]. With /(n, . . . , r,-) 
being given by the l.h.s. of Eq. (|102J) . the sum on the 
r.h.s. of Eq. I|10()|l contains classes of (2 — o~2j)j terms 
that give the same result (note that for j = 2 the cyclic 
permutation and its reverse coincide, so we have only j 
instead of 2j terms in the class). By forming a set P(j) £ 
P(j) containing exactly one representative of each class 
(where obviously P(j) has j!/ [(2 — S2j)j] members), the 
sum can thus be simplified, leading to Eq. 1(55)1. 

With this preparation at hand, we may derive Eq. 
JSU [together with Eqs. JS^J) and Ipfy] by following these 



steps: We substitute Eqs. gJJ and lO into Eq. l(T3|). 
multiply out and perform all spatial integrals over delta 
functions, where only terms of the form S(sj — s^+i) con- 
tribute [terms of the form 5(ta— Sj) giving zero integrals, 
because of n(r A ) = 0, cf. the remark below Eq. (|47|l ] . 
After renaming the remaining integration variables ac- 
cording to 



/(si, ...,s/) 



n [ d3s i n ( s i 
=1 e fn/d'Vfc) 



(i)' 



M0> 



(103) 



7reP(0 b=i • 
the result may be written in the form 

k 



A k U A {v A ) = Y. A k U A(rA) 



(104) 



1=1 



with 



1 



A l k U A (r A )=- J ,\„ 



II/: ■- 



d 3 Sj n(s 



x ^2 Q Vl (si,iu)---q ni (si,iu) 

J)i>0,...,t7i>0 
»;i + ...+rji=fe — ( 

xU AB ...B(rA,Si,...,Si,iu), (105) 



where each power of the factor 



|e 1 n(r)a B (w) 



1 - ±e 1 n(r)a B (u) 



(106) 



is due to the integration of one term containing 5(sj — 
Sj+i), and 

f 00 

/ duU AB ... B {rx,...,r l+1 ,iu) 
Jo 

= U AB ... B {ri,...,ri +1 ), (107) 

recall Eq. (|52p. Summing Eq. I|104|) over k, and rearrang- 
ing the double sum, we find 

00 00 k oo 

£ A k U A (r A ) =J2Y1 A i U ^A) = £ A l U A (r A ), 

(108) 



k=l 

where 



k=l 1=1 



1=1 



-1 />OC * /> 

A l U A (r A ) = - du / d 3 Si 



n(sj) 



1 



±£ 'llfSjjQjjH %=Q 

x C/AB...B(ryi,si, . . . ,si,iu). 



q^fa,™) 



(109) 
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After performing the geometric sums 

oo 

Xy(r,w) - 1 - \e^n{v)a B {uj), (110) 
j=o 

cf. Eq. i(106l) . the denominators in Eq. (|109fl cancel, so by 
recalling Eqs. ifTTj). (|107|1 and JTUHJ, we arrive at Eq. (|5TJ 
together with Eqs. ((SI and (|5^|l . 



B Calculation of the two-point correlation term 
for the dielectric ring [Eqs. ()66jl and (|67[)] 

An approximation to the two-point correlation term in 
the retarded limit as given by Eq. l(41(l - ((44[) [together 
with Eqs. dHHI— (ESI) and Eq. in the case of the di- 

electric ring can be obtained by replacing the variable 
si by its average across the cross section of the ring 
(\ta — si | ~ pa for a -C r ), evaluating the si-integral, 
and separating the S2-integral into two parts, 

a) _ ^(0)>'(0) 

{p\a pa p2-rr p27T—Xa/ro j 

/ dz / dp p / d<f> + ira 2 / r d6> > 

J -\a JO JO Jxa/r j 

1 



For the integral in A^Ua^a), the approximations 
a = 1 ~p A , /3~2r o |sm(0/2)|, 
a $ = /3-7 : 



r o |sin(0/2)| 



7-q: 



PA 

2r 2 sin 2 (0/2) 
Pa 



1 



(116) 



are valid for a <C ro . Inspection of Eq. 1111(1 shows that 
the leading term in (a/ro) of zA^C/a^a) is due to the 
factor /3 3 oc sin 3 (6>/2) in the denominator of the integrand 
[cf. Eq. H118|) below], and comes from regions where 
sin(0/2) <C 1. Hence we may apply a Taylor expansion in 
powers of sin(0/2), retaining only 



/i(a,/?,7) 



13 



/ 2 (/3,7,a) 



3 
i' 



(117) 



Substituting Eqs. I|116|) and l|117|) into Eq. (TTTJ), and 

performing the 0-integral using 



2tt — Aa/ro 



d(9 



/ Aa/ro |sin^/2)| 
eventually leads to 



8 



■o[ln(Aa/r )], (118) 



X a^(a+p +1 ) [McM+Mw P) (d • /3) 2 

+ / 2 (a,/3, 7 )(^-7) 2 + /2(/3, 7 ,a)(7-a) 2 
+ / 3 (a,A7)(a-/3)(M)(7-a)] 
= Z^' c t/ A (r A ) + ^[^(r^) (111) 

where the integral in A^' c Ua{ya) extends over an ap- 
proximately cylindrical volume of cross section ira 2 and 



2 7hcVa A (0)x 2 (0) v 1 

5127r J eoPyi ^ 



so that 



7/Wa A (0)x 2 (0) 



5127r 3 £ PA 



x /(A), 



where 



/(A) 



2A 



1 

A 2 ""' VTTa 2 " 



(119) 



(120) 



(121) 



volume of the remaining open ring. 



length 2Aa, and that in A 2 Ua{?a) extends over the [recall Eq. (|115|l ]. Note that the approximations made 

for calculating A 2 ,c Ua(ya) break down for large A while 
those made for calculating A 2 ' t Ua(ya) break down for 
small A. We put 



For the integral in A 2 ' c Ua(ta), we may approximate 

a = 7 ~ p A , (3 ~ \/ z 2 + p 2 , 



c,$ = -0.^J^M=, 7 .d.-l (112) 

yj + p Z 

for a<r , and Eqs. d2)-(HJ [recall Eq. lpT7|)] simplify 
to 

13 15 

fi(a,f3,j) ~ y, / 2 (a,/3,7) = f 2 (i,a,(3) ~ y, 

MAT,")--, / 3 (a,/3, 7 )^-y. (113) 

Substituting Eqs. (|TT2^ and (fTT5| into Eq. (jTTTl) . carry- 
ing out the ^-integral, and using 



/(A) 



± 7 



max /"(A) + min f (A) 

0.5<A<1.5 0.5<A<1.5 



max /(A) — min /(A) 

0.5<A<1.5 0.5<A<1.5 



3.5 ± 1.4, 



(122) 



in Eq. 1120(1 . resulting in Eq. ftffijl. 

A similar procedure may be applied in the nonre- 
tarded limit, where Eq. 1(46(1 leads to 



A 2 U A {r A ) = 



3HV 



Xa , r , 2z 2 - p 2 2A 
dz / dpp — = — ^=^= , 

-Xa Jo y/z 2 + p 2 V 1 + A 2 



one may find 



A 2 2 c U A (r A ) = 



7hcVa A (0) X 2 (0) A 



256^e p 7 A s/T+J 2 ' 



(114) 



(115) 



128^%o Jo 

rXa pa />2n 

< / dz I dpp / d(j) + ira 2 



'-Xa JO JO 

x l-3(d-/3)(/3- 7 )( 7 .d) 

a 3 /3 3 7 3 

^' c tf A (r A ) + A 2 /U A (r A ) 



ductA(iu)x (iu) 

2-rr — Aa/ro 
Aa/ro 



r o d0 



(123) 
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Use of Eqs. I|112|) and l|114|) leads to 
3W 



A 2 2 > c U A (p A ) = 



64tt 3 £oP^ Jo 
A 

x , 



duaA(iu)x (iu) 



The spatial integrals in Eqs. (JHSJ [recall Eq. (|%4l ]. 
(|87|l and H88|) can be carried out explicitly for p(z) = 
cos 2 (k z z), 



(124) 



while using Eq. (|116[) . neglecting the term (a ■ 0)((3 
7) (7 • a), and recalling Eq. I|118|) . results in 



d "^» 82 «->-3ra)> 

/•OO />oo 

/ dzcos 2 (k z z) dz' e~ 2bz ' cos 2 (k z z') 
2b 6 + 8b 4 k 2 + 5b 2 kf + 2fc6 



A^Ua(pa) = 



dua A (iu)x (iu) x 



86 2 (& 2 + fc 2 ) 2 (6 2 + 4/c 2 ) . 



(131) 



(132) 



A 2 

(125) 

Combining Eqs. (|124f) and l|125|l in accordance with 
Eq. 1|123[> . we obtain 

AjZ/aGcu) = 777^ r / dua A {iu)x 2 (iu) x /(A) , 



resulting in 



A\U A (z A ) = ~TS~ 9 / duu 2 a j4 (m)x(m) / d<? 

2 , / x 2 



128^ 3 e ^ 



x e 



-262A 



26 2 



(126) 



b 2 + k 2 



bcY Ifu 
~u)~ + 2\b~c 



b 

(133) 



which, in combination with Eq. 023, implies Eq. fTJ. A\U A (z A ) = [°°du u 2 a A (iu) X 2 (m) [ dq 1 

16tt - 2 J Jo b 



C Asymptotic power laws in the case of a half 
space with oscillating susceptibility [Eqs. (|9()[) 
and (I9B1 



As a preparing step, we derive the linear and quadratic 
expansions in x of the coefficients 



C 4 = 



3hca A (0) 
647r 2 e 

' 2 I 



d „, [ [x(Q) + ifr - VxWT^ 
1 w \ [ x (o) + 1]« + vx(o) + 7; 



g - a/x(0) + v 2 x J_ 

V + y/x{0) + V 2 W 4 



and 



A 1 C 4 + A 2 C 4 



(127) 



167r 2 e J 

:A 1 C 3 +A 2 C 3 



dw a J 4(zu)— — 



X(*u) + 2 



(128) 



x e 



-26z A 



26 2 



A 2 U A (z A ) 



b 2 + k 2 







3 








V U J 


"4 H 







x e 



-262A 



duu ceA(iu)x (iu) / dg- 
16tt 2 7 y b 

2b 6 + 8b 4 k 2 + hb 2 k\ + 2k 6 z 



, (134) 
(I 



(6 2 4 


VI 


2 (6 2 + 


4fc 2 ) 




1 






2~{te)~ 


F 4 







(135) 



In analogy to the procedure outlined in Ref. [221 , the 
retarded limit may conveniently be treated by introduc- 
ing the new integration variable v — bc/u, transforming 
integrals according to 



duu 2 I dq?-e- 2bzA 



c 3 f°°dv r dbb3e -2 bZA 

Jl Jo 



(136) 



that can be found for the retarded and nonretarded dis- 
tance laws of the homogeneous half space [231 - Substi- 
tuting 



applying the approximation l|23(l . and carrying out the 
u-integrals. Application of this procedure to Eqs. (|133|) . 
(TgUl , and l(T35j) leads to 



[X(0) + l]v - y/x(0) + v 2 
[ X (0) + l]v + ^/x(0)+v 2 



1 1 

2 ~ iv 2 



v-^ X (0)+v 2 



x(0)~ 



1 



1 1 

4 ~ 8^ 



x 2 (o) 



1 



« + Vx(0) + v 2 4w 



= -x3x(0) + 3-rx 2 (0) 



8v 4 



(129) 
(130) 



into Eq. (|127|l and carrying out the remaining w-integral, 
we arrive at Eqs. and JjJIjJl, while Eq. (|128|l implies 
Eqs. Eg ) and ^5t . 



4„ j( , ) = ^M«!wW, (138 , 



6407T 2 e 4 



A 2 U A {z A ) 



A3hca A (0)x 2 {0)H 3 (k z z A ) 
89607T 2 e 4 



(139) 



where we have introduced the definitions l)96|) and H97|) . 
Combining Eqs. p37|l - (|139|l in accordance with Eq. (JJSJl 
and using Eqs. (J22J) and ((^3), we arrive at Eq. J5U). 
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The asymptotic behavior of Eqs. 1)133}! 1)135)1 in the 

nonretarded limit may be obtained by transforming the 
integral according to 

/>oo />oo 

/ du dq\ e- 2bZA . . . 
Jo Jo b 

rOO r OO 

i > / du dbe- 2bZA ... , (140) 
Jo Jo 

retaining only the leading power of u/(bc), carrying out 
the ^-integral and discarding higher-order terms in uza/c 
(cf. Ref. ), resulting in 

A\U A (z A ) = - H ^ zZA J [°° dua A (iu) X (iu), (141) 

A\U A {z A ) = H ^ zZ H I™ dua A {iu) X 2 {iu), (142) 

A 2 TT i \ hnoH (k z z A ) f°° 2 2 
A 2 U A (z A ) = — — / duu a A (iu)x {m) , 

6ZTT Z A J 

(143) 

recall Eqs. flUEl and ((HZJ) - Upon using Eq. f7R|. Eqs. (JHU 
- 1)143)1 lead to Eq. 191)1 . where we have neglected the 
term proportional to zT 1 in consistency with the nonre- 
tarded limit and used Eqs. <|M)l and (|9*3)l . 
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